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ABSTRACT 

We present self-consistent cosmological magnetohydrodynamic (MHD) sim- 
ulations that simultaneously follow the formation of a galaxy cluster and the 
magnetic field ejection by an active galactic nucleus (AGN). We find that the 
magnetic fields ejected by the AGNs, though initially distributed in relatively 
Ck small volumes, can be transported throughout the cluster and be further ampli- 

fied by the intra-cluster medium (ICM) turbulence during the cluster formation 
process. The ICM turbulence is shown to be generated and sustained by the 
frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process 
is shown to exist in the ICM and amplify the magnetic field energy and flux. The 
total magnetic energy in the cluster can reach ~ 10 61 ergs while micro Gauss (/iG) 
fields can distribute over ~ Mpc scales throughout the whole cluster. This finding 
shows that magnetic fields from AGNs, being further amplified by the ICM tur- 
bulence through small-scale dynamo processes, can be the origin of cluster-wide 
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1. Introduction 

There is growing evidence that the ICM is permeated with magnetic fields, as indicated 
by the detectio n of large-scale, diffused radio emission c alled radio halos and relics (see re- 



cent reviews by lFerrari et al.ll2008l ; ICarilli fc Taylorll2002l ). The radio emissions are extended 



over > 1 Mpc, covering the whole cluster. By assuming that the total energy in relativistic 
electrons is comparable to the magnetic energy, one often deduces that the magnetic fields in 
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the cl uster halos ca n reach 0.1 — 1.0 fiG and the total magnetic energy can be as high as 10 61 
ergs (IFerettil Il999l ). The Faraday rotation measurement (FRM), combined with the ICM 
density measurements, often yields cluster magnetic fields of a few to ten /j,G level (mostly 
in the cluster core region). More interestingly, it reveals tha t magnetic fields can have a 
Kolmogorov-like turbulent spectrum in the cores of clusters JVoet & Enfilinl l2003h with a 
peak at several kpc. Other studies have suggested t hat the coherence scales of magnetic 
fields can range from a fe w kpc to a few hundred kpc (jEilek fc Owen! 120021 ; iTaylor fc Perley 
19931 ; IColgate fc Lill2000l ). implying large amounts of magnetic energy and fluxes. Under- 
standing the origin and effects of magnetic fields in clusters is important because they play a 
crucial role in determining the structure of clusters through processes such as heat transport, 
which conse quently aff ect the applicability of clusters as sensitive probes for cosmological 
parameters (IVoitll2005l ). 



Although the existence of cluster-wide magnetic fields is clear, their origin is still poorly 
understood. Two scenarios have received most a ttention: a) magnetic fields are initially 
from the outflows of norm al or active galaxies (IDonnert et al.l 120091 ; iFurlanetto fc Loeb 
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Kronberg et al. 



20011 ). and such fields can be further amplified by cluste r merger 



jRoettiger et alill999h and turbulence JDolag et al.ll2002l : bubois k Tevssier!l2008l); and b) 
very s mall proto-galact ic seed fields are amplified by dynamo processes in clusters ( IKulsrud et al 
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Rvu et al.l 120 08). though the exact mechanism for dynamo is still being debated 



(IBernet et al.ll2008l ). Large scale radio jets from AGNs serve as one of the most intriguing 



candidate s in the first scenario because they could carry large amount of magnetic energy 
and flux (jBurbidgd Il959l ; iKronberg et al.l l200ll ; ICroston et al.l 120051 ; iMcNamara &: Nulsen 
20071 ). The magnetization of the ICM and the wider inter-galactic medium (IGM) by AGNs 



has been suggested o n the energetic grounds ([Colgate fc Lil |2000| ; IFurlanetto fc Loebl 12001 



Kronberg et al.ll200ll ). though the exact physical processes of how magnetic fields will be 
transported and amplified remain sketchy. Some of the key physics questions in these mod- 
els include: what is the volume filling process of AGN magnetic fields in the ICM/IGM? 
Is the ICM turbulent and what are its properties? Is there a dynamo in the ICM that can 
amplify fields? 

In this Letter, we present self-consistent cosmological MHD simulations to address the 
question of the origin of magnetic fields in clusters. We explore specifically the scenario 
that cluster-wide magnetic fields initially came from the magnetic fields of an AGN. We also 
describe the properties of the ICM turbulence and demonstrate the existence of dynamo in 
the ICM. The simulations are described in Section [2] and results are presented in Section [31 
Discussions are given in Section 4. 
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Basic Model and Simulations 



We have performed detailed cosmological MHD simulations of galaxy cluster formation 
with the magnetic field injection from an AGN, using the newly developed ENZO+MHD 
code, which is a n Eulerian cosmological MHD code with adaptive mesh refinement (AMR) 
( IXu et al.ll2008l ). Our simulations use a ACDM model with parameters h = 0.7, Q m = 0.3, 
Qb = 0.026, = 0.7, and a$ = 0.928, with initial conditions extracted from the Simulated 
Cluster Archive (see http://lca.ucsd.edu/data/sca/). The simulation volume is 366 Mpc on 
a side, and it uses a 128 3 root grid and 2 level nested static grids in the Lagrangian region 
where the cluster forms. This gives an effective root grid resolution of 512 3 cells (0.5 Mpc) 
and dark matter particles of mass 1.49 x 1O 1O M . 

The simulations were evolved from redshift z = 30 to z = adiabatically. We "turn on" 
the AGN magnetic injection at redshift z — 3, centered at a massive galaxy in a proto-cluster 
which has a virial radius r v ~ 0.15 Mpc, a baryon virial mass « 3.2 x 10 11 M , and 
a virial total mass m v ~ 7.6 x 10 12 M & . The cluster eventually grows to r v ~ 2.15 Mpc, 
rrib pa 9.0 x 10 13 M , and m v ~ 1.1 x 10 15 M by z = 0. The total injected magnetic energy 
by AGN is about 2 x 10 60 ergs, with an average input power of 1.75 x 10 45 ergs s -1 for a 
duration of 36 Myr. The magnetic energy is injected inside 0.2 r v . Because it is currently 
not possible to resolve both the galaxy cluster and the AGN environment simultaneously, we 
have adopted a n approach that mimics the possible magnetic energy injection by an AGN 
( ILi et al.ll2006l ). The size of the injection region and the associated field strength are not 
realistic when compare d to the real AGN jet s, but on global scales, the previous studies by 
Nakamura et al.l (120061 ) and Ku et al.l (120081) showed that this ap proach can reproduce the 
observed X-ray bubbles and shock fronts ( McNamara et al. 20051 ). AMR is allowed only in 
a region of (50 Mpc) 3 where the galaxy cluster forms. During the cluster formation, the 
refinement is controlled by baryon and dark matter overdensity. In addition, all the region 
where magnetic field strength is higher than 10 -7 G will be refined to the highest level. There 
are a total of 8 levels of refinement beyond the root grid, for a maximum spatial resolution 
of 11.2 kpc. Consequently, our simulations are equivalent to ~ 600 3 uniform grid MHD runs 
in the cluster region with full cosmology. The simulation was performed on the linux cluster 
Coyote at LANL with about 300,000 CPU hours consumed. 
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3. Results 

3.1. Global Morphology 

To illustrate the formation and evolution of the cluster along with the evolution of 
magnetic fields from the AGN, we present images of the projected gas density and magnetic 
energy density at different stages of cluster formation in Fig. [TJ At z = 3 (t = 0), an AGN 
is "turned on" in a sub-cluster (which eventually merges with another sub-cluster about 
400 Myr after the AGN injection has finished). The AGN is "turned off" at t = 36 Myr 
(z = 2.95). At t — 180 Myr, we see the density cavitie s produced by t he AGN magnetic fields, 



reminiscent of the jet-lobe structure of radio galaxies (IXu et al.ll2008l ). At a later time t = 468 
Myr, the jet-lobe structure is destroyed by the cluster mergers. When two sub-clusters finish 
merging at about t = 1.1 Gyr, a large part of magnetic fields is carried out of the cluster 
center region by the ejected flow from mergers. As the evolution progresses, magnetic fields, 
which follow the plasma motion, are being sheared, twisted, and spread throughout the whole 
cluster. Judging by the images from t — 1.1 to t — 6.28 Gyr, this volume-filling process is 
quite efficient (see Discussions). By z = 0.5 (t = 6.28 Gyr), magnetic fields are well mixed 
with the ICM and are distributed throughout the whole cluster, with some high magnetic 
field regions obviously from shock compressions. At z = 0.0, while the cluster has relaxed, 
the magnetic fields seem to distribute over the whole cluster quasi-uniformly. 



3.2. Energy Evolution and Magnetic Field Radial Profile 

The evolution of the total thermal, kinetic, and magnetic energy inside the cluster's 
virial volume is shown in the top panel of Fig. [2J The kinetic energy is calculated as the 
turbulent component by subtracting the bulk flow motion. By the AGN injection, ~ 2 x 10 60 
ergs of magnetic energy is input into a sub-cluster. A significa nt fraction (~ 8 0%) of this 



energy is deposited into the ICM due to expansion and heating ( jXu et al.l 120081 ). When the 
two big sub-clusters merge at t — 400 Myr, the total thermal and kinetic energies of the 
cluster increase by ten-fold by t ~ 600 Myr. During this major merger, the magnetic fields 
are still very local, largely in magnetic bubbles. So the merger event itself did not significantly 
change the magnetic field energy. The thermal and kinetic energies of the cluster continue to 
increase as the cluster grows in mass via accretion and mergers as indicated by the variations 
in the energy evolution curves. At t ~ 2 Gyr, the magnetic fields from the AGN have been 
spread throughout a significant volume of the cluster (see Fig. [T]). Starting from this time, 
the magnetic energy experiences an exponential increase by a factor of 20 until t ~ 6 Gyr, 
at which time the growth has slowed. After t ~ 6 Gyr, the cluster grows slowly in its total 
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energy and becomes relaxed, and the magnetic energy increases slowly with the growth of 
the cluster. By z — 0, the total magnetic energy has reached ~ 10 61 ergs inside the cluster. 
Considering that the magnetic energy will drop as oc r~ l due to the the expansion of the 
universe, the total magnetic energy has actually increased by ~ 75 times between t = 2 and 
t = 10 Gyr. 

In the bottom panel of Fig. [2j we present the spherically averaged radial profiles of 
magnetic field strength at different epochs. At z = 0, the magnetic field strength is ~ 1.5 
fiG at the core, decreasing slowly to ~ 0.7 fiG at ~ Mpc radius. These radial profiles are 
different from wha t is shown by other simulations which often ex hibit faster radial decline 
(IDolag et al.l 120021 ; iDubois &; Teyssierl l2008t iDonnert et al.l 120091 ) . It is presently unclear 
whether this difference is caused by the different origins of "seed" magnetic fields or by the 
effects of different numerical techniques and resolutions. This difference deserves further 
study since it may be used to distinguish different origins of magnetic fields in clusters. 
Furthermore, it will be nec essary to compare our simulation result s with cluster magnetic 
observations in detail (e.g. 
in future publications. 



Govoni et al.ll2006l ; iGuidetti et al.l 120081 ). This will be presented 



3.3. Small-scale Turbulent Dynamo and MHD Turbulence 

The physical origin for the exponential amplification of the magnetic energy is due to 
a cluster-wide turbulent dynamo process. In Fig. [3] we present the power spectra of the 
ICM plasma's kinetic energy density and magnetic energy density in a comoving cube with 
5.71 Mpc on the side. The scales are shown in the comoving units so that they are not 
affected by the universe expansion and the power changes are from the dynamics of cluster 
evolution alone. Here, k = 0.003 kpc -1 corresponds to the maximum radius of the central 
region 2.86 Mpc where a high spatial refinement is adopted. Since the cluster's total thermal 
energy is a factor of ~ 3 — 5 larger than its kinetic energy, the ICM can be thought as nearly 
incompressible plasmas overall, though shocks generated by mergers are very frequent and 
important. In fact, the flows and shocks from mergers t end to be global and propagate 



through the whole cluster (see also the earlier work by iRoettiger et al.l 119991 ) . We can 
divide the full kinetic spectrum approximately into three regions: the "injection" region for 
k ~ 0.003 — 0.01 kpc -1 where the large scale flows and shocks from the mergers "stir" the 
whole cluster; the "cascade" region for k ~ 0.01 — 0.1 where the spectrum shows a smooth 
power law similar to a Kolmogorov-like incompressible turbulence; the "dissipation" region 
for k > 0.1 where the spectrum steepens gradually. We attribute this steepening to both 
the dissipation by shocks (which has a physical origin) and the limited spatial resolution 
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(which has a numerical origin). These features in the kinetic energy density spectrum lead 
us to conclude that the ICM is turbulent and our simulations have captured the essence of 
this turbulence. This ICM turbulence is in a driven-dissipative state where frequent mergers 
will drive the turbulence over relatively short time scales (a few Myr) but the turbulence is 
decaying in-between mergers on timescales of ~ Gyr (which is approximately the dynamic 
timescale for the whole cluster). 

The magnetic energy density spectrum also has three regions in the A;— space corre- 
sponding to the kinetic energy spectrum. From t — 2 to t = 10 Gyr (z = 1.5 to z = 0.1), the 
magnetic spectrum for k ~ 0.01 — 0.6 retains an "invariable" shape but the energy density 
undergoes exponential amplification then goes i nto saturation. This is a strong sig nature 
for the so-called small-scale turbulent dynamo ([Brandenburg & Subramanianl 120051 ). Fur- 
thermore, the magnetic energy density peaks at k ~ 0.2 with ~ 3 x 10 -16 ergs cm" 3 . The 
corresponding kinetic energy density is ~ 8 x 10~ 16 ergs cm -3 . So, the magnetic energy for 
k > 0.1 has saturated at a sub-equipartition level (by a factor of ~ 3). The drop-off at high k 
should be due to the numerical dissipation. For k < 0.01, neither kinetic or magnetic energy 
seems to have saturated. It is interesting to note that the dynamo process starts vigorously 
only at t ~ 2 Gyr (z = 1.5), when the magnetic fields have been spread through a significant 
fraction of the whole cluster (see Fig. [Q. Putting Figs. [1] - [3] together, we see that the 
ICM turbulence both amplifies the magnetic energy and diffuses the fields throughout the 
cluster. The magnetic energy density saturation occurred at t ~ 6 Gyr but magnetic fields 
continue to "grow" in their spatial extent through turbulent diffusion. Furthermore, the 
un-signed magnetic flux through the mid-plane of the cluster has exponentially grown from 
7.6 x 10 41 to 1.23 x 10 43 G cm 2 from z = 2 to z = 0.5, which is another clear indication of 
the turbulent dynamo that is responsible for both amplifying the field energy and diffusing 
the field through the cluster. 



3.4. Faraday Rotation Measurement 



We have also computed the synthetic FRM by integrating to the mid-plane of the cluster. 
Fig. H] shows the spatial distribution of FRM at z = 0. The typical value of FRM is ±200 rad 
m~ 2 , with high values concentrated in the cluster core region. Interestingly, the FRM map 
not only shows the small scale variations reminiscent of the ICM MHD turbulence, but also 
displays long, narrow filaments with dimensions of 300 kpc x 90 kpc. The FRM magnitudes 
and spatial distribu tions from simulations are quite consistent with observations o f radio 
galaxies in clusters iGuidetti et alil2008l : lEilek fc Owenll2002l ; iTavlor fc Pe"rlevlll993h . 
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Discussions 



The results presented here need to be taken as an initial step in better understanding 
the evolution of AGN magnetic fields over the lifetime of a cluster. The fact that the AGN 
magnetic fields are injected very early in the history of the cluster formation (i.e., z=3), 
before the major merger event s, could be imp ortant. For AGN fields injected relatively late 
in the clus ter formation ( e.g., Ku et al.l 120081 1 or injected into relatively quite background 
ICM (e.g.. lLiu et al.ll2008l ). the magnetic fields may not experience extensive stretching and 
transport so that they could survive in a bubble morphology without undergoing significant 
mixing with the background ICM. 

Note that even though we have solved the ideal MHD equations, there is clearly nu- 
merical diffusion that has allowed the magnetic fields to diffuse in the ICM. The rate of 
diffusion is often related to the numerical Reynolds number and magnetic Reynolds number. 
We estimated that these numbers in our simulations are on the order of a few hundred (this 
is especially true in the cluster core region where most of the magnetic field en ergy resides). 
The small scale dynamo theory and simulations (jBoldyrev fc Cattaneol 120041 ) have shown 
that the dynamo will grow under such conditions, w hich is consistent with our findings. The 
real ICM could have Reynolds number as low as ten ([Reynolds et al.ll2005l ). but its magnetic 
Reynolds number is less well determined, especially in a magnetized turbulent medium. 

Since AGNs are commonly observed in galaxy clusters, one important implication of 
our studies is that the magnetic fields from AGNs alone are perhaps enough to "seed" the 
cluster and the ICM turbulence generated and maintained by mergers will amply and spread 
the AGN fields via a dynamo process. This is consistent with observa tions that clu sters with 
large scale radio emissions are often correlated with cluster mergers (IFerettil 120051 ). 
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Fig. 1. — Snap shots of the projected baryon density (upper rows) and magnetic energy 
density (lower rows) for different epochs of cluster formation and evolution. Each image 
covers a region of 5.71 Mpc x 5.71 Mpc (comoving). The projected results are obtained 
by integrating 5.71 Mpc (comoving) centered at the cluster along lines perpendicular to the 
observed plane. The eight columns are marked with the time t since the AGN injection and 
the respective redshift z. The top panel uses different color scale for each plot to have the 
best visual effect. The color range of the bottom panel is the same for all subplots, ranging 
from 4.2 x 10 19 to 1.3 x 10 22 cm -2 for the baryon particle number density and from to 5 
x 10 11 ergs cm -2 for the integrated magnetic energy density. 
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Fig. 2. — Top panel: Temporal evolution of different components of energy inside the virial 
radius of the cluster. The variations in the thermal and kinetic energies are due to mergers. 
Bottom panel: The spherically averaged radial profile of magnetic field strength at different 
epochs of the cluster formation. The radius is measured in the proper frame. They show the 
turbulent diffusion of magnetic fields throughout the cluster, yet maintaining its strength 
via the dynamo process. 
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Fig. 3. — Power spectra of the kinetic energy density and magnetic energy density of the 
ICM at different epochs. The ICM turbulence is represented by the Kolmogorov-like spectra 
in kinetic energy. The magnetic energy is amplified via a dynamo process. 
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Fig. 4. — Faraday rotation measurement of the ICM by integrating to the mid-plane of the 
cluster. It covers a region of 2.86 Mpc x 2.86 Mpc (comoving) at z=0. The color range 
shown is from —200 (blue) to 200 (red) rad m -2 . The peak value of rotation measurement 
is about ±400 rad m~ 2 . Filamentary structures are quite common. 



